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Abstract 

The widths of a glueball decay to two pious or kaons are analyzed in the 
pQCD framework. Our results show that the glueball ground state has small 
branching ratio for two-meson decay mode, which is around 10~^. The pre- 
dicted values are consistent with the data of ^ — > vrvr, KK if ^ particle is a 
2++ glueball. Applicability of pQCD to the glueball decay and comparison 
with XcJ decay are also discussed. 
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I. INTRODUCTION 



The existence of glueballs and hybrids is a direct consequence of QCD. Since the J/ip 
is predicted to have appreciable decays to '-ygg where two gluon formation is expected to 
enhance the production of tensor and scalar glueballs, radiative J/ip and ip' decays have 
long been realized to be a favored area for glueball searches. At present one pays particular 
attention to three states: /o(1500), /j(1710) (J=0 or 2), and ^(2230) (J > 2). They are 
glueball candidates of O"*"*" or 2++ states. 

In particular, the BES collaboration discovered new, non-strange decay modes of the 
^(2230) state, such as ^ ^ tttt and ^ J^P HI- Compared with other mesons, the ^(2230) 
has many distinctive properties [0: 

(1) Flavor-symmetric decays to light hadrons. After removal of the phase space factor, 
the probability for ^ — »• 7r"'"7r~ is of the same order as that for ^ — >■ K^K~ . 

(2) Copiously production in radiative J /ip decays. From the upper limit of 1 x 10~^ for 
Br{^ —>■ pp)Br{C, —>■ KK) P,^, where KK includes all kaon pairs, a lower bound 3 x 10^^ 
for the Br{J/il) — > 7,^) can be estimated. 

(3) Narrow width. Both results from Mark III and BES show that the width of the 
^(2230) is only about 20 MeV. Assuming = 20 MeV, one can easily estimate from (2) 
that the Br{^ — > K^K~) and Br{!^ — >■ tt+tt^) are smaller than 2%, resulting in partial 
widths F^+^- and V^+k- smaller than 400 KeV 0. 

As a consequence, the qq model, multi-quark model and hybrid model can not easily 
reproduce three observations above. On the other hand, these properties are naturally 
explained by identificating ^(2230) as a glueball state with J^^ = 2++. 

Since the observed decay modes into vrvr, KK and pp are expected to be only a small 
portion of the decay modes of the ^(2230), searches for other decay modes are very impor- 
tant. From theoretical point of view, the narrow width of a glueball follows from a loose 
interpretation of the OZI rule. The gluons in the glueball would annihilate and create a qq 
pair that would form the lighter hadrons. Since this suppression only acts at one vertex, it 
is so called vOZI rule For example, the total width of a 2~^~^ glueball can be estimated 
from VOZI rule to be about ^^^^(mo)!^"^ — 20 MeV. 

In order to understand the narrow decay width quantitatively, we study a pure 0++ or 
2++ glueball decay to two light mesons in perturbative QCD (pQCD). Our results show that 
a pure glueball decay to two light mesons has a small branching ratio or a narrow width, 
and this conclusion can be generalized to be valid for any two mesons. As a consequence, 
there is no dominant decay channel for a pure 0+^ or 2"*"+ glueball. 

The paper is organized as follows. As a comparison with the glueball, the formula for 
Xcj decay to two mesons in pQCD is reviewed briefly in Sec. II. Glueball decays into two 
mesons are calculated in Sec. III. In Sec. IV the numerical results and applicability of pQCD 
to glueball decays are discussed. The last section is reserved for summary and conclusions. 

II. HEAVY QUARKONIUM WITH J^^ = 0++ OR 2++ 

It is interesting to compare a glueball decay with the Xcj- A brief review for Xcj decay 
will be given in the pQCD framework. For a heavy quarkonium coupling to gluons, the 
vertex can be obtained with its radial wave function at origin, or its differentiation. Since 







the wave function is sharply peaked at small internal momentum for heavy quarkonium, the 
nonrelativistic limit is a good approximation. 

As an example, we review the derivation of Bt{xcJ tctc) in pQCD. To leading order, 
the Fermi movement of quarks in pion can be neglected comparing with the momentum flow 
of order m^, the mass of c quark, in the hard scattering amplitude. Thus, the amplitude of 
Xcj decay to two pions can be expressed in a factorizated form as vr form factor as 

A = J dxdy(f)^{x)THix, y, Q^)(l)^{y), (1) 

where 4>n{x) is the distribution amplitude of pion obtained by integrating the transverse mo- 
mentum of the Bethe-Salpeter wave function 0. The hard scattering amplitude Th{x, y, Q^) 
(see fig. 1) will include information of a heavy quarkonium coupling to two gluons: 
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with li and I2 the momentum of gluons. 
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and the wave function of heavy quarkonium in nonrelativistic quark model (NRQM) 
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X{p,k) = Y.{27r)6{k'' - —){LM;SS,\JJ,)^LM{k)PssMk), 
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where ipLAiik) and PssXVi ^) spatial and spin part, respectively. T^^ involve the spin part 
of pion wave function and present the coupling of gg — > tttt: 



(5) 



For quarkonium in P-wave state, only 0"''^ and 2++ can decay to two pions. In the case 
of 2++ state, the polarization tensor e""^ is composed of the spin and orbital polarization 

vectors, e'^{Sz) and e'^(M), of Xc2 as 



(6) 



The polarization tensor satisfies e"^ = e'^", p^e"'^ = 0, e'^a = and 



where 
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Now we arrive at the decay width if as{xiX2Q'^) — 
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where = ^ is color factor and the integral I2 is given by 



dx(f)^{x) f dy(t)^{y) 
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Similarly, The 0^^ decay width can be expressed as 
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with the overlap integral Jq given by 



dx(t>n{x) f dy(t>^{y) 
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The total width of Xcj can be obtained by calculating the width of its decay to two real 
gluons in a similar way. In NRQM, it is 
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Finally, we get the branching ratio 
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which is independent of Xcj wave function. 

Numerical calculations can be easily done. However, the obtained results are much 
smaller than the data Even when stretching all parameters to their extreme values, 
the predictions stay a factor 3-6 below the data [|10|. Moreover, the obtained branching 
ratio dominantly come from the end-point region, where the pQCD is not available. As 
stated by Bodwin, Braaten, and Lepage |Tl|, the color-octet decay contribution arising 
from the higher Fock component \ccg) of the Xcj wave function is actually not suppressed 
with respect to the usual color-singlet decay, owing to the angular momentum barrier in 
|cc) contribution, and is necessary to separate rigorously short-distance effects and long- 
distance effects. Including the color-octet contributions, branching ratio of Xcj tttt can 
be compared with the experimental data |]T0 |. 
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III. GLUEBALL 



Similar to the QQ state, we apply pQCD to the glueball decay. The amplitude for 
its decay to two pions can be written as eq. (|I|), too. Now contain the glueball wave 
function. The coupling of a state to two gluons can be obtained from the requirements 
such as linear in the polarization vectors of two gluons, ei and £2, and Lorentz and gauge 
invariant P,p!2[. Keeping only the leading twist term, the wave function can be written 
phenomenologically as 
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where ^ is color wave function, e^'^ is the polarization tensor of 2++ glueball and V^'^ is 

1,2), e* and kl is the polarization vector and 
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defined in eq. 
momentum of the i-th gluon, respectively. 

The leading order contribution to a 2"*"+ glueball decay to two pions in pQCD is shown 
in fig. 2. In center of mass frame, the amplitude can be written as 

A = glj dzM^)dz'Mz')Mk)TL E eM.(^)G;,(sl)G^p(^2)e^(sl)e;(52)^'^^ (17) 

J ^1^2 S1S2 

where T^^ is the same as in eq. (|^) and the color factor will be included in the width 
formula. The momenta of gluons are fixed by that of quarks in pions, i.e., ki = zQi + z'Q2 
and k2 = 'zQi + 'z'Q2- The width can be obtained 



^(2^ 



TTTT 



c 



TT 



30M 



dzdz'(j)Tr{z)4>n{z')<f)2{ki)as — - 



zz + z'z' 



1^1 



ZZZ'Z 



where C is the color factor | and z = [1 — z). The total width of 2^ 
its decay width to two on-shell gluons with k'j = 
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can be obtained by 
(19) 



Similarly, we can obtain the width of 0^+ decay to two pions 
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and the total width 
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Comparing eqs. (|T8|) and (pO]) with eqs. (|9|)-([T3|), it can be shown that there are two 
different ingredients between the glueball decay and the P-wave quarkonium decay, i) There 



is no end-point singularity in the expression (|T^) and (^0]) for the glueball decay since gluons 
in the glueball are directly related to the glueball wave function; ii) Since we get the total 
width by the decay to two on-shell gluons, it depends on the glueball wave function at a 
particular momentum point. Especially, the point is located at k"^ = Intuitively, 
the wave function must be peaked at low momentum, since a composite particle has little 
amplitude for existing while its constituents are flying apart with large momentum. While 
the mass M is large, the obtained width decreases fast and depends on the shape of the 
wave function drastically. On the other side, the next order contributions to the total width, 
whose amplitudes are expected to have the form as J d^k(j){k)f{k) in general, will be lack of 
the suppression of the wave function and not sensitive to the shape of the wave function as 
long as f{k) is smooth. It means that the higher order contributions are important to the 
total width if the mass is large, since the zero-th order contribution is suppressed strongly 
by wave function. In another word, the pQCD evolution of the wave function to a particular 
scale will be nontrivial. 

The formula for decay to a pair of kaons can be easily obtained by substituting the 
distrubtion amplitude of pions by that of kaons. As a result of isospin symmetry, the 
branching ratio to ttOttO is half of the charged channel except that 71%^ in the final state 
can be formed via QCD anomaly. This contribution is negligible due to the minor difference 
between the mass of u-quark and cZ-quark. The anomaly contributions to rjr], rjr]' and rj'r]' 
channels are not necessarily small, since two soft gluons have large possibility to form a rj 



[t]') meson |jT3 



{^\^-^Gl^G'^nri) = f^f.K (22) 

However, it is a low energy theorem. If the energy scale of the decay process is high enough, 
the anomaly contribution is also negligible as it emerges via QCD renormalization and is 
essentially a higher order contribution in Ug- For example, the branching ratio of Xcj VV 
has no apparent enhancement comparing with vr^vr^ channel. So the branching ratios of 7]7], 
rji]' and rj'T]' channels are expected to have the same order of tttt channel if pQCD dominate 
the decay process. 



IV. APPLICABILITY OF PQCD AND NUMERICAL RESULTS 



Evaluation of eqs. (|T8|~^T]) will require the glueball wave function, which we have little 
information about up to now. In the following context we will use an oscillator wave function 



(f){k) = Agexp(-6jA;^ 



(23) 



with parameter appropriately chosen. If 0^"*" and 2^~^ glueballs have similar wave functions 



a constraint on bg can be obtained from their total 



in momentum space as argued in ref. |T4 
decay width. The 0+"*" glueball candidates /o(1500) and /o(1700) have widths 100 ~ 150 
MeV. The 2^~^ glueball candidate ^(2230) has width ~ 20 MeV. From our width formulae 
(|g) and (^ we have 



^(2230) _ 4M|A2 exp(-262M|) 20 MeV 



ry„(i5oo) l5MlAlexp{-2blMl) 120 MeV 



(24) 



If the masses of scalar and tensor states are equal, the ratio will be simply ^, which is in 
accord with ref. |jl2|]. Taking into account the mass gap, we get 6^ ~ 0.4 GeV~^ if A2 — A^. 



In general, Ag is larger for a higher mass state. For example, Ak is larger that (see 
below). So we get bg > 0.4. According to the discussion in the previous section, the higher 
order corrections may be large and should be larger for a higher mass state. Thus we can 
say that the lower limit of 6^ is 0.4 GeV~^. Furthermore, the parameter bg is related to the 
radius of hadron. We know that b^ ~ 0.8 GeV~^ for a well established pion oscillator wave 
function |15]. Assuming that all hadrons have similar size, 6^ will be in the same range for 
glueball. The parameter 6^ for the glueball will be taken in the range 0.6 < 6^ < 0.8 in the 
following context. 

Before getting the numerical results it is necessary to discuss whether the derived formula 
is applicable or reliable. The applicability of the pQCD theory to exclusive processes at the 
present experimental energy region was argued by Isgur and Smith [|16| by excluding the 
contributions of the end-point regions where sub-leading (higher twist) terms are a priori 
likely to be greater than the perturbative contribution. Recently, the applicability of pQCD 
to the pionic form factor has been examined by cutting end-point contribution and 



by including the effects of Sudakov form factor [|T^. The first approach argues that the 



pQCD results are self-consistent in the energy region where the contributions after cutting 
dominate. The second one attempts to explain the suppression in the end-point regions by 
including the effects of Sudakov form factor of the quarks, which serves as a natural filter 
to pick out the hard contributions. Two approaches give similar conclusion that pQCD is 
applicable to the pion form factor as > 4 GeV^. In fact, the applicability of pQCD 
to exclusive processes differs from one process to another and it depends on the end-point 
singularities. For example, the hard scattering amplitude of the XcJ decay to two pions 
is more singular than the case of the pion form factor. The obtained rate of Xcj tttt 



depends strongly on the end-point behavior of the pion wave function |T9|. In the case 
of the glueball decay, the hard scattering amplitude of glueball decay to two pions is less 
singular than that of the pion form factor, and it has a good behavior to ensure that the 
dominant contributions come from the hard part. We can argue that pQCD is applicable for 
the glueball, particularly, for the glueball candidate ^(2230). We will adopt the first method 
which is simpler. 

In order to realize the condition as{zzM'^) < 1 in eqs. (|T8|) and (pOD, we extend the 
parameterization of as{Q^) by replacing pO[| 

127r 

""^^^'^ ^ (33-2n;)ln(g2 + M2)/A2 ^^^^ 

to reflect the fact that at low the transverse momentum intrinsic to the bound state 
wave-function flow through all the propagators. The parameter Mq can be determined to 
ensure as{0) = 1. The cut contribution, obtained by cutting the integral limit to keep only 
hard contributions, say, < 0.7, will be compared with the uncut contribution. 

In order to calculate the decay width we take the wave function of the pion and the kaon 

as 
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where ntq means the mass of m or d quark and uis the mass of s quark. The distribution 
amphtude (t){x) can be obtained by integrating fc^: 
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The parameters can be adjusted by using the constraints from decays of vr — 77 and 
TT — > /ii/ {K — ^ iip) and the average quark transverse momentum {k'])^, ~ {k\)K — (360 
MeV)2 [|1|,|21[. It turns out 



mq = 0.25 GeV, 



0.80 GeV"^ A^ = 27.7 GeV" 



0.55 GeV, = 0.80 GeV'^ Ak = 51.0 GeV" 



The numerical results are listed in the table 1. The parameter bg is taken as 0.6, 0.7 and 
0.8. Here we haven't fixed the normalization constant and the branching ratio is independent 
of Al. The branching ratios of 2++ and 0"*"+ decay to tt+tt" or K^K~ are around 10~^ for 
a pure glueball. The numerical result is consistent with ^ — > tttt if ^ is a 2"'""'" glueball. 

Also we examine the applicability of pQCD to the glueball decay by cutting end-point 
contributions. Our numerical results show that the hard part (a^ < 0.7) contribute about 



11% to ^(2^ 



TT ' TT 



), whose value is listed in table 1, and the ratio is nearly the same for 



different 6^ because the partial width is not sensitive to 6^. For 0+"*" the ratio is about 50%. 
Another kind of distribution amplitude of pion, the asymptotic form, is also used. We find 
that the branching ratio increases about 10% for 2+"*" glueball and 20% for 0^"*" one and the 
hard contributions occupy about 61% in the width of 2++ decay to tt+tt^ and 30% in the 
width of 0+"*" decay. 



V. CONCLUSIONS AND SUMMARY 



We have analyzed the decay width of a glueball decay to two pions and kaons in the 
pQCD framework and it can be generalized to other two-meson decay channels. The nu- 
merical results of the decays show that the branching ratio is small and the decay width is 
very narrow, compared to the qq bound state. The branching ratios are consistent with the 
data of ^ — 27r, 2K if ^ particle is a 2"*"+ glueball. 

Our conclusions follows: 

(1) Any glueball candidate produced from J/ip radiative decay has small branching 
ratio for two-meson decay mode, which is around 10~^. This result is derived in the pQCD 
framework. However, if the mass of glueball is not high enough to ensure that pQCD 



s 



dominate the decay process, the QCD anomaly will play an important role and the 7777, rjr]', 
and rj'ri' channels will enhance apparently. 

(2) Applicability of pQCD to the glueball decay is discussed in our paper. We show that 
the hard scattering amplitude of glueball decay has a good behavior at the end-point region 
and is favored by pQCD while XcJ is opposite. The reason is that gluons are directly related 
to the wave function for glueball decay. The extra quark propagator (and the differentiation 
on it) produces the end-point singularity in Xcj exclusive decay. More specific speaking, 
pQCD is applicable for ^(2230) decay to two mesons if it is a 2++ glueball, but is not 
reliable for /o(1500). 

(3) Related to the end-point singularity in XcJ decay, the color octet contributions have 
the same order as color singlet and may be the essential part for P-wave state quarkonium 
decay. Therefore, the dynamic mechanism will be very different for decays of glueball and 
quarkonium with the same J^^. 

(4) It is not difficult to generalize our calculation to other two-meson decay processes. 
The conclusion is expected to be similar. Therefore the branching ratio is small for each 
decay mode and there is no dominant decay channel for a pure glueball whose mass is high 
enough to ensure the pQCD contributions dominate. 
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FIGURES 

FIG. 1. The hard scattering diagram for XcJ ~^ t^t^- 
FIG. 2. The hard scattering diagram for the gluebah decay to tttt. 



TABLES 



TABLE L Numerical results. 







^tot 




i?r(7r+7r 




Br{K+K~) 




0.6 


2.82x10-^^2 


1.08x10-^^2 


0.38x10- 


-2 


0.69x10-2 


^(2230) 


0.7 


1.04x10-^^2 


1.01x10-'^A2 


0.97x10" 


-2 


1.75x10-2 




0.8 


3.86x10-*'A2 


0.95x10-'^A2 


2.45x10" 


-2 


4.45x10-2 




0.6 


8.46x10^^^2 


6.00x10-6^2 


0.71x10" 


-2 


1.11x10-2 


/o(1500) 


0.7 


5.39x10-^^2 


5.64x10-6^2 


1.05x10" 


-2 


1.65x10-2 




0.8 


3.44x10-4^2 


5.32x10-6^2 


1.55x10" 


-2 


2.47x10-2 



19 




Fig. 1 




Fig. 2 



